Multi-agent systems, Equiprobability, Gamma distributions and 

other Geometrical questions 



Ricardo Lopez-Ruiz 
DIIS and BIFI, Facultad de Ciencias, 
Universidad de Zaragoza, E-50009 Zaragoza, Spain 

Jaime SanudoQ 

Departamento de Fisica, Facultad de Ciencias, 
Universidad de Extremadura, E-06071 Badajoz, Spain 

Xavier CalbeJl] 
Instituto de Astrofisica de Canarias, 
Via Ldctea, s/n, E-38200 La Laguna, Tenerife, Spain 
(Dated: February 1, 2008) 

Abstract 

A set of many identical interacting agents obeying a global additive constraint is considered. 
Under the hypothesis of equiprobability in the high-dimensional volume delimited in phase space 
by the constraint, the statistical behavior of a generic agent over the ensemble is worked out. The 
asymptotic distribution of that statistical behavior is derived from geometrical arguments. This 
distribution is related with the Gamma distributions found in several multi-agent economy models. 
The parallelism with all these systems is established. Also, as a collateral result, a formula for the 
volume of high-dimensional symmetrical bodies is proposed. 

PACS numbers: 87.23. Ge, 05.90.+m, 89.90.+n 
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I. INTRODUCTION 



In this paper, a general multi-agent system under an additive constraint is considered. 
Following the derivations done in Refs. Illj2|j3|, we work out the geometrical properties of 
this system in phase space in order to obtain its statistical behavior. We observe a striking 
coincidence. The different dynamical mechanisms that have been proposed in the literature 
to model the interaction among agents in multi-agent economic systems^^^ display the 
same statistical results than those derived from the geometrical properties of our system. 
This fact seems to suggest a close relationship between the local interactions among the 
agents in the former systems and the global geometrical conformation in phase space of our 
general system. We are inclined to think that all those dynamical mechanisms provoke that 
those systems, with an adequate change of coordinates, evolve equiprobably over the volume 
of accessible states in the transformed phase space. 

We start in Section II by recalling the particular results obtained in Ref. 3J for the cases 
in which the constraint has a linear or quadratic dependence on the variables. Then, in 
Section III, we obtain the statistical behavior for a more general constraint. In Section IV, we 
establish a possible relationship with other systems^^ 1 ^ in which the Gamma distributions 
are also obtained. A formula for the volume of high- dimensional symmetrical bodies is 
proposed in Section V. The last Section VI contains our conclusions. 



II. RECALLING SOME RESULTS 

(A) Let us assume N agents interacting in an open economy^, each one with coordinate 
Xi, i = 1, . . . , N, with Xi > representing the wealth or money of the agent i, and a total 
available amount of money E. The additive constraint reads: 

Xl+X2^ h £jv_i + X N < E. (1) 

Under random evolution rules for the exchanging of money among agents^, let us suppose 
that this system evolves in the interior of the A-dimensional pyramid given by Eq. ([T]). 
We can suppose that the state or the bank system of western societies plays in this model 
the role of a heat reservoir that supplies money instead of energy. The formula for the 
volume Vn{E) of an equilateral A-dimensional pyramid formed by N + 1 vertices linked by 
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iV perpendicular sides of length E is 



E N 

V N (E) = — . (2) 



If each point on the iV-dimensional pyramid is equiprobable, then the probability f(xi)dxi 
of finding the agent % with money Xi, with normalization condition f(xi)dxi = 1, is 
proportional to the volume formed by all the points into the (N — l)-dimensional pyramid 
having the ith-coordinate equal to Xj. We have shown 3 that f(xi) verifies 

If we call e the mean wealth per agent, E = Ne, then in the limit of large N (N — ► oo), we 
have 

^ VN t { f~ X) =-e^dx, (4) 
n»i V N (E) e K ' 

where the index % has been removed because the distribution is the same for each agent, and 

thus the wealth distribution can be obtained by averaging over all the agents, 

f(x) = e" 1 e~ x/e dx. (5) 

This Boltzmann-Gibbs distribution has been found to fit the real distribution of incomes in 
western societies^. 

(B) Now let us suppose a one-dimensional ideal gas of iV non-identical classical particles 
with masses rrii, with i = 1,...,N, and total maximum energy E. If particle i has a 
momentum rriiV j, we define a kinetic energy: 

K = p 2 t = -niiVi, (6) 

where pi is the square root of the kinetic energy Then the constraint reads: 

P 2 i+pI + ---+P 2 n-i+P 2 n<E. (7) 

We see that the system has accessible states with different energy, which is supposed to 
be supplied by a heat reservoir. These states are all those enclosed into the volume of the 
iV-sphere given by Eq. ([7j), with radius E 1 / 2 . The formula for the volume Vn(R) of an 
iV-sphere of radius R is 

JV 

Vn(R) = fTJ^Tj^' ( § ) 



where T(-) is the Gamma function. If we suppose that each point into the iV-sphere is 
equiprobable, then the probability f(pi)dpi of finding the particle i, with coordinate pi 
(energy pf) and normalization condition J^ R f(pi)dpi = 1, is proportional to the volume 
formed by all the points into the (N — l)-sphere having the zth-coordinate equal to pi. We 
have shown^ that f{pi) verifies 

nP '> Vn(E^) ' [ ' 

If we call e the mean energy per particle, E = Ne, then in the limit of large N (N — > oo), 
we have 

where the index i has been removed because the distribution is the same for each particle. 
Thus the asymptotic distribution 

/(p) = /J^ 1/2 e- p2/2e (11) 



can be obtained by averaging over all the particles. If the change of variables p = v is 
performed, with v the generic velocity of a particle, then the Maxwellian distribution is just 
derived from geometrical arguments. 



III. MULTI-AGENT SYSTEMS AND EQUIPROBABILITY: 
GENERAL DERIVATION OF THE ASYMPTOTIC DISTRIBUTION 

In this section, we address the same problem above presented but in a general way. Let 
b be a positive real constant (cases b = 1,2 have been indicated in the former section). If 
we have a set of positive variables (xi, x 2 , ■ ■ ■ , x^) verifying 

xl + xl + '-' + x^ + x^^E (12) 

with an adequate mechanism assuring the equiprobability of all the possible states 
(xx, X2, ■ ■ ■ , xn) into the volume given by expression (|T2l) . will we have for the generic variable 
x the distribution 

f(x)dx ~ e- l l h e~ xbtht dx, (13) 

when we average over the ensemble in the limit N, E — ► oo, with E = Ne, and constant e?. 
Now we show that the answer is affirmative. 
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From the cases b = 1,2, (see Eqs. ([3]) and (jHJ)), we can extrapolate the general formula 
that will give us the statistical behavior f(x) of the generic variable x, when the system runs 
equiprobably into the volume defined by a constraint of type (1121) . The probability f(x)dx of 
finding an agent with generic coordinate x is proportional to the volume Vn-i({E — x h ) l l b ) 
formed by all the points into the (N — l)-dimensional symmetrical body limited by the 
constraint (E — x b ). Thus, the iV-dimensional volume can be written as 

V^E 1 ^) = f El/ " V N ^((E - x b ) l ' b ) dx. (14) 
Jo 

Taking into account the normalization condition 1 f(x)dx = l, the expression for f(x) 
is obtained: 

fix) - Y^ii^m (15 ) 

The ^-dimensional volume, Vjv(6, p), of a 5-symmetrical body with side of length p is 
proportional to the term p N and to a coefficient gb{N) that depends on iV: 

V N (b,p)=g b (N)p N . (16) 

The parameter b indicates the original equation (fl2l) that defines the boundaries of the 
volume Vjv(6, p). Thus, for instance, from Eq. fl2]), we have gb=i(N) = 1/N\. 

Coming back to Eq. ffTol) . we can manipulate Vn((E — x b Y^ b ) to obtain (the index b is 
omitted in the formule of Vjv): 

N_ 

r i N « / T \ b 

V N {{E - x b Y' b ) = g b {N) [{E-x b yi b ] = g b (N) E* 1 1 - . (17) 

If we suppose E = Ne, then e represents the mean value of x b in the collectivity, that is, 
e =< x b >. UN tends toward infinity, it results: 



limfl-^l = e~ xb l b \ (18) 

iV»i \ E ) y 1 

Thus, 

V N ((E - x b ) l/b ) = V N (E 1/b ) e~ xb/be . (19) 

Substituting this last expression in formula f|T5|) . the exact form for f(x) is found in the 
thermodynamic limit [N,E — > oo): 

f{x)dx = c b e- 1/b e- x "/ be dx, (20) 



b\ h 



with Cb given by 



c b ~ gb{N ~ l) (21) 



Hence, the conjecture ( 1131) is proved. 

Doing a thermodinamical simile, we can calculate the dependence of e on the temperature 
by differentiating the entropy with respect to the energy. The entropy can be written 
as S = — kNf£° f(x) In f(x) dx, where f(x) is given by Eq. (120]) and k is the Boltzmann 
constant. If we recall that e = E/N, we obtain 

. kN f E\ kN , 
S(E) = — In + — (1 - blncb), (22) 

where it has been used that e =< x b >= J °° x b f(x)dx. 
The calculation of the temperature T gives 

m , fdS\ kN k 

r i4=B=s' (23) 

Thus e = kT/b, a result that recovers the theorem of equipartition of energy for the quadratic 
case 6 = 2. The distribution for all b is finally obtained: 

1/6 



f(x)dx = c b (J^) e x " /kT dx. 



(24) 



IV. GAMMA DISTRIBUTIONS 



If we perform the change of variables y = e x l h x in the normalization condition of f(x), 
Io° f( x )dx = 1, we find that 

" poo . "I — 1 

c b = / c-y' b dy . (25) 



If we introduce the new variable z = y /b, the distribution f(x) as function of z reads: 

f(z)dz = -^r z 1 ^ 1 e~ z dz. (26) 

b b 

Let us observe that the Gamma function appears in the normalization condition, 

f '«*=Af ^'^iMD" 1 - (27) 

This implies that 

,-^f y (28, 



6 



By using Mathematica the positive constant c& is plotted versus 6 in Fig. 1. We see that 
linife^o Q> = oo, and that lim^oo c& = 1. The minimum of c& is reached for 6 = 3.1605, 
taking the value c& = 0.7762. Still further, we can calculate from Eq. (l2"g|) the asymptotic 
dependence of Cf, on b: 

limc^^^(l-l + ...), (29) 

lim c b = b- l ' b f 1 + 2 + . . .V (30) 

where 7 is the Euler constant, 7 = 0.5772. The asymptotic function (l29l) is obtained 
after substituting in (|28|) the value of T(l/6) by (1/6—1)!, and performing the Stirling 
approximation on this last expression, knowing that 1/6 — > 00. The function (IHUj) is found 
after looking for the first Taylor expansion terms of the Gamma function around the origin 
x = 0. They can be derived from the Euler's reflection formula, r(x)r(l — x) = ir / sin(7ra;). 
We obtain T(x — > 0) = + T'(l) + • • •. From here, recalling that T'(l) = —7, we get 
T(l/6) = 6 — 7 + • • •, when 6 — > 00. Although this last term of the Taylor expansion, —7, 
is negligible we maintain it in expression (1301) . The only minimum of q, is reached for the 
solution 6 = 3.1605 of the equation ip(l/b) + log6 + 6—1 = 0, where ?/>(■) is the digamma 
function (see Fig. 1). 

Let us now recall two interesting statistical economic models that display a statistical 
behavior given by distributions of the form (126]) . that is, the standard Gamma distributions 
with shape parameter 1/6, 

f(z)dz= f Lz*- 1 e-'dz. (31) 



MODEL A: The first one is the saving propensity model introduced by Chakraborti 
and Chakrabarti^. In this model a set of N economic agents, having each agent % (with 
i — 1, 2, • • • , N) an amount of money, Ui, exchanges it under random binary (i, j) interactions, 
(ui,Uj) — > (u'ijUj), by the following the exchange rule: 

u'i = Xui + e{l - X)(ui + Uj), (32) 
u'j = Xuj + e(l — X)(ui + Uj), (33) 
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with e = (1 — e), and e a random number in the interval (0,1). The parameter A, with 
< A < 1, is fixed, and represents the fraction of money saved before carrying out the 
transaction. Let us observe that money is conserved, i.e., U{ + Uj = u\ + u'j, hence in this 
model the economy is closed. Defining the parameter n(A) as 

n(X) = y^y, (34) 

and scaling the wealth of the agents as z = nuj < u >, with < u > representing the average 
money over the ensemble of agents, it is found that the asymptotic wealth distribution in 
this system obeys the standard Gamma distribution 6 

f(z)dz = zzz — r z n - 1 e~~ z dz. (35) 

The case n = 1, which means a null saving propensity, A = 0, recovers the model of 
Dragulescu and Yakovenko^ in which the Gibbs distribution is observed. If we compare 
Eqs. (1351) and (}3~Tj) . a close relationship between this economic model and the geometrical 
problem solved in the former section can be established. It is enough to make 

n = 1/b, (36) 
z = z, (37) 

to have two equivalent systems. This means that, from Eq. fl36|) . we can calculate b from 
the saving parameter A with the formula 

b = I— * (38) 
1 + 2A v ; 

As A takes its values in the interval (0, 1), then the parameter b also runs in the same interval 
(0,1). On the other hand, recalling that z = x b /be, we can get the equivalent variable x 
from Eq. fl33), 



x 



-a/ 



u 



<u> 

,b 



(39) 



where e is a free parameter that determines the mean value of x in the equivalent geometrical 
system. Formula fl39l) means to perform the change of variables Ui — > Xi, with i — 1, 2, • • • , N, 
for all the particles/agents of the ensemble. Then, we conjecture that the economic system 
represented by the generic pair (A, u), when it is transformed in the geometrical system given 
by the generic pair (b, x), as indicated by the rules (13"B1) and (I3"9~j) . runs in an equiprobable form 
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on the surface defined by the relationship (jT2l . where the inequality has been transformed 
in equality. This last detail is due to the fact the economic system is closed, and then it 
conserves the total money, whose equivalent quantity in the geometrical problem is E. If 
the economic system were open, with an upper limit in the wealth, then the transformed 
system would evolve in an equiprobable way over the volume defined by the inequality (1121) . 
To see more clearly the equivalence between surface and volume in a statistical ensemble we 
address the reader to Ref. [jj where this question has been discussed for similar cases. 

MODEL B: The second one is the model introduced by Angle^. In this model a set 
of N economic agents, having each agent i (with i = 1, 2, • • • , N) an amount of money, u^, 
exchanges it under random binary interactions, (ui,Uj) — > (u'^u'j), by the following the 
exchange rule: 

Ui = Ui - Am, (40) 
u'j = Uj + Au, (41) 

where 

Am = r](xi — Xj) euxi — [1 — r](xi — Xj)] euxj, (42) 

with e a random number in the interval (0, 1). The exchange parameter, u, represents the 
maximum fraction of wealth lost by one of the two interacting agents (0 < u < 1). Whether 
the agent who is going to loose part of the money is the z-th or the j-th agent, depends 
nonlinearly on (xj — Xj), and this is decided by the random dichotomous function r](t): 
r]{t > 0) = 1 (with additional probability 1/2) and r](t < 0) = (with additional probability 
1/2). Hence, when x% > Xj, the value rj = 1 produces a wealth transfer from agent % to agent 
j with probability 1/2, and when Xi < Xj, the value rj = produces a wealth transfer from 
agent j to agent % with probability 1/2. Defining in this case the parameter n{u) as 

. . 3 — 2ui , 
»M = — , (43) 

and scaling the wealth of the agents as z = nuj < u >, with < u > representing the average 
money over the ensemble of agents, it is found that the asymptotic wealth distribution in 
this system obeys the standard Gamma distribution 8 

f(z)dz = —— z n - x e" 2 " dz. (44) 
1 \n) 



The case n = 1, which means an exchange parameter u = 3/4, recovers the model of 
Dragulescu and Yakovenko^ in which the Gibbs distribution is observed. If we compare 
Eqs. (1441) and (131 p . a close relationship between this economic model and the geometrical 
problem solved in the last section can be established. It is enough to make 

n = 1/6, (45) 
z = z, (46) 



to have two equivalent systems. This means that, from Eq. (145]) . we can calculate b from 
the exchange parameter uj with the formula 

b = (47) 

As uj takes its values in the interval (0, 1), then the parameter b runs in the interval (0, 2). 
It is curious to observe that in this model the interval w £ (3/4, 1) maps on b G (1, 2), a fact 
that does not occur in MODEL A. On the other hand, recalling that z = x b /be, we can get 
the equivalent variable x from Eq. (|46|) . 

1/6 



X 



(48) 



.< u > 

where e is a free parameter that determines the mean value of x b in the equivalent geometrical 
system. Formula (j48p means to perform the change of variables Ui — > Xi, with i — 1, 2, • • • , N, 
for all the particles/agents of the ensemble. Then, we conjecture that the economic system 
represented by the generic pair (A, u), when it is transformed in the geometrical system given 
by the generic pair (b, x), as indicated by the rules (1471) and (148j) . runs in an equiprobable form 
on the surface defined by the relationship (fT2"l) . where the inequality has been transformed 
in equality. As explained above, this last detail is due to the fact the economic system is 
closed, and then it conserves the total money, whose equivalent quantity in the geometrical 
problem is E. If the economic system were open, with an upper limit in the wealth, then 
the transformed system would evolve in an equiprobable way over the volume defined by the 
inequality ( fl2i) . This equivalence in using the surface or the volume of a statistical ensemble 
in order to obtain its statistical behavior has been discussed in Ref. [j for similar cases. 
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V. OTHER GEOMETRICAL QUESTIONS 



We shall proceed now to derive an asymptotic formula (N — > oo) for the volume of the 
N- dimensional symmetrical body enclosed by the surface 

x\ + x\ + • • • + x h N _ x +x b N = E. (49) 

The linear dimension p of this volume, i.e., the length of one of its sides verifies p ~ E x l h . 
As argued in Eq. f|T6|) . the A-dimensional volume, VnQ), p), is proportional to the term p N 
and to a coefficient <7&(iV) that depends on A. Thus, 

V N (b,p)=g b (N)p N , (50) 



where the characteristic b indicates the particular boundary given by equation ff49|) . 
For instance, from Eq. ([2]), we can write in a formal way: 



N 



9>=i(N) = r( | + 1) . (51) 

From Eq. ([8]), if we take the diameter, p = 2R, as the linear dimension of the iV-sphere, we 
obtain: 

^ F{f+T)' (52) 

These expressions ( l5Ti) and ( l52i) suggest a possible general formula for the factor gb{N), let 
us say 

JV 

»w - Fj^ry (53) 

where a is a 6-dependent constant to be determined. For example, a = 1 for 6=1 and 
a = 7r/4 for 6 = 2. 

In order to find the dependence of a on the parameter b, the regime N — > oo is supposed. 
Applying Stirling approximation for the factorial (y)! in the denominator of expression (|53j) . 
and inserting it in expression ( 12TI) . it is straightforward to find out the relationship: 

c b = (ab)- 1 ^. (54) 

From here and formula ff28l) . we get: 

r /1 \i & 

(55) 



r[i + i 



n 



that recovers the exact results for b = 1,2. The behavior of a is monotonous decreasing 
when b is varied from b = 0, where a diverges as a ~ 1/6+ ■ • •, up to the limit b —>■ oo, where 
a decays asymptotically toward the value a m = e~ 7 = 0.5614. 
Hence, the formula for gb{N) is obtained: 



It would be also possible to multiply this last expression (|56l) by a general polynomial K(N) 
in the variable N, and all the derivation done from Eq. ( |53i) would continue to be correct. 
We omit this possibility in our calculations. For a fixed N, we have that gb{N) increases 
monotonously from gb{N) = 0, for b = 0, up to <?&(-/V) = 1, in the limit b — > oo (see Fig. 2). 
For a fixed b, we have that gb{N) decreases monotonously from gb{N) = 1, for iV = 1, up to 
gb(N) = 0, in the limit iV — > oo (see Fig. 3). 

The final result for the volume of an iV-dimensional symmetrical body of characteristic b 
given by the boundary (1491) reads: 



with p ~ E x l h . 

VI. CONCLUSIONS 

In this work, we have considered a general multi-agent open system verifying an addi- 
tive constraint. Its statistical behavior has been derived from geometrical arguments. The 
Maxwellian and the Boltzmann-Gibbs distributions are particular cases of this type of sys- 
tems. Also, other multi-agent economy models, such as the Chakraborti and Chakrabarti's 
mode£, the Angle's model^, and the Dragalescu and Yakovenko's model^, show similar sta- 
tistical behaviors than our general system. This fact suggests a geometrical interpretation 
of all those models. The equivalence with the Chakraborti and Chakrabarti's model is es- 
tablished when the geometrical characteristic b of our model runs in the interval (0, 1). The 
equivalence with the Angle's model is established when b varies in the interval (0,2). As a 
particular case of both types of model, the Dragulescu and Yakovenko's model is obtained 
for 6 = 1. 




(56) 




(57) 
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Let us remark that we have not found in the literature other multi-agent models to 
establish an equivalence with our system in the range b G (2, oo). This point remains as 
an open question and it can be a challenge that will probably trigger other works in this 
direction. 
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Figures 
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FIG. 1: Normalization constant q, versus b, calculated from Eq. (j28|) . The asymptotic behavior 
is: limft—jo Q> = and lim^^oo c& = 1. This last asymptote is represented by the dotted line. The 
minimum of q, is reached for b = 3.1605, taking the value q, = 0.7762. 
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FIG. 2: The factor gb{N) versus b for N = 10,40,100, calculated from Eq. (|56|) . Observe that 
g b (N) = for b = 0, and lim 6 ^oo 9b{N) = 1. 




FIG. 3: The factor g b (N) versus N for b = 10,40,100, calculated from Eq. (J55]). Observe that 
5 6 (AT) = 1 for N = 1, and lim iV -*ooS&(A r ) = 0. 
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